Measurements of cosmic-ray energy spectra with the 2""^ CREAM flight 
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During its second Antarctic flight, the CREAM (Cosmic Ray Energetics And Mass) baUoon experiment cob 
lected data for 28 days, measuring the charge and the energy of cosmic rays (CR) with a redundant system of 
particle identiflcation and an imaging thin ionization calorimeter. Prehminary direct measurements of the abso- 
lute intensities of individual CR nuclei are reported in the elemental range from carbon to iron at very high energy. 



1. Introduction 

The CREAM experiment was designed to mea- 
sure the composition and energy spectra of cos- 
mic rays approaching energies up to lO'^^ eV. 
Since December 2004, four instruments were suc- 
cessfully flown on balloons over Antarctica where 
they collected several million CR events in the el- 
emental range from hydrogen to iron, and with 
total particle energies reaching the 100 TeV scale 
and above. The final goal of the experiment is to 
provide a deeper understanding of the accelera- 
tion mechanism of cosmic rays and to test the va- 
lidity of the astrophysical models describing their 
propagation in the Galaxy [1] . 
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In this paper, we present preliminary energy spec- 
tra of the even-charged, abundant nuclei from car- 
bon to iron as measured by the instrument during 
its second flight (CREAM-II). 

2. The CREAM-II instrument 

The instrument for the second flight included: 
a redundant system for particle identification, 
consisting (from top to bottom) of a timing- 
charge detector (TCD), a Cherenkov detector 
(CD), a pixelated silicon charge detector (SCD), 
and a sampling imaging calorimeter (CAL) de- 
signed to provide a measurement of the energy of 
primary nuclei in the multi-TeV region. 
The TCD is comprised of two planes (120 x 
120 cm^) of four 5 mm-thick plastic scintillator 
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paddles, read out by fast timing photomultiplier 
tubes (PMT). The CD is a 1 cm-thick plastic ra- 
diator, with 1 surface area, read out by eight 
PMTs via wavelength shifting bars. The SCD is 
a dual layer made of 312 silicon sensors, each seg- 
mented as an array of 4x4 pixels, which covers 
an effective area of 0.52 m^ with no dead regions. 
The CAL is a stack of 20 tungsten plates (50x50 
cm^ , each 1 Xq thick) with interleaved active lay- 
ers of 1 cm-wide scintillating fiber ribbons, read 
out by 40 hybrid photodiodes (HPD). It is pre- 
ceded by a 0.47 Ai„t-thick graphite target to in- 
duce hadronic interactions of CR nuclei. 
The second CREAM payload was launched on 
December 16*'' 2005 from McMurdo and flew over 
Antarctica until January 13*'' 2006, at a balloon 
float altitude between 35 and 40 km. A detailed 
description of the instrument and its flight per- 
formance can be found in 

3. Data analysis 

A data set collected in the period Dec. 19*''- 
Jan. 12*'', under stable instrument conditions, 
was used in this analysis. The first step of the 
procedure consists of selecting events with an ac- 
curate trajectory reconstruction of the primary 
particle, and measuring its charge and energy. 
The direction of the particle is given by the axis of 
the shower reconstructed in the imaging calorime- 
ter. This is obtained by a fit of the candidate 
track points sampled along the longitudinal de- 
velopment of the shower. On each CAL plane, 
they are defined as the center of gravity of the 
cluster formed by the cell with maximum signal 
and its two neighbours. The fitted shower axis is 
back-projected to the target and its intersections 
with the two SCD layers define the impact points 
of the primary particle, with a spatial resolution 
better than 1 cm rms. In each SCD plane, the 
pixel with the highest signal is located inside a 
circle of confusion with a 3 cm radius centered on 
the impact point. 

If the signals of the matched hits in the upper 
and lower SCD planes are consistent within 30%, 
they are selected as two independent samples of 
the specific ionization dE/dx and corrected for 
the pathlength traversed by the particle in the 



300 






250 






C 

> 200 




I 






Number of 










so 


\ 

t 1 




i 

1 1 1 





6 8 10 12 14 16 18 20 22 24 26 28 

Charge (Z) 



Figure 1. Charge histogram obtained by the SCD 
in the elemental range from boron to iron. 

silicon sensors. A charge Z is assigned to the 
particle by averaging the two measured values of 
the specific ionization and taking into account its 
dependence. The charge distribution recon- 
structed by the SCD is shown in figure [U by fit- 
ting each peak to a gaussian, a charge resolution 
cr is estimated (in units of the electron charge e) 
as: 0.2 for C, N, O; ~ 0.23 for Ne, Mg, Si; ~ 0.5 
for Fe. A 2cr cut around the mean charge value is 
applied to select samples of each element, while 
for iron a la cut is used. In this way, 1288 O are 
retained, as well as 456 Si and 409 Fe candidates. 

The total energy [Ed), deposited in the 
calorimeter by an interacting nucleus, is measured 
by summing up the calibrated signals of its cells. 
In order to infer the primary particle energy E 
from Ed^ an unfolding procedure is applied. In 
fact, due to the finite energy resolution of the 
detector, the measured counts in a given energy 
interval must be corrected for overlap with the 
neighbouring bins. This can be done by inverting 
the matrix equation: Mi = ^ijNj where Nj 
and Mi are the "true" and measured counts in 
each energy bin, respectively. A generic element 
of the mixing matrix Aij represents the probabil- 
ity that a CR particle, carrying an energy cor- 
responding to a given energy bin j, produces an 
energy deposit in the calorimeter falling in the 
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Figure 2. Energy deposited in the calorimeter 
by a selected sample of carbon nuclei. Simulated 
(histogram) and real (dots) events are shown. 



bin i. A detailed MonteCarlo (MC) simulation 
of the instrument, based on the FLUKA 2006.3b 
package |5 , was developed to estimate the unfold- 
ing matrix. Sets of nuclei, generated isotropically 
and with energies chosen according to a power- 
law spectrum, are analyzed with the same proce- 
dure as that used for the flight data. Each ma- 
trix element A^j is calculated by correlating the 
generated spectrum with the distribution of the 
deposited energy in the calorimeter [J]. In order 
to get a reliable set of values of the unfolding ma- 
trix for each nucleus, the MC simulation is finely 
tuned to reproduce both flight data and the cal- 
ibration data collected with accelerated particle 
beams [5] . The agreement of the MC description 
with the real instrument behaviour was carefully 
checked. As an example, in figure [5] the response 
of the calorimeter to carbon nuclei from the fiight 
data is compared with an equivalent set of simu- 
lated events. 

4. Energy spectrum 

The interval of energy spanned by the mea- 
sured CR events is divided into bins of width AEi, 
larger than the energy resolution of the calorime- 
ter and centered at values E™''^''' (calculated ac- 
cording to the definition of [6]). For each element, 
the "true" number of counts Ni in each bin is 
obtained as a result of the unfolding algorithm. 
The absolute differential intensity $ at an energy 
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AEi e X TOI X TOA x ST^ x T 



where T is the exposure time, e the efficiency of 
the selection cuts, Sft the geometric factor of the 
instrument, and TOI and TOA are, respectively, 
the corrections to the top of instrument and to 
the top of the atmosphere. 

The geometric factor Sil is estimated from MC 
simulations by counting the fraction of generated 
particles entering the upper SCD plane and cross- 
ing the upper CAL plane. A value of 0.46 m^sr 
is found. During the flight, the live time T was 
measured by the housekeeping system onboard. 
The selected set of data amounts to a live time 
of 16 days and 19 hours, close to 75% of the real 
time of data taking. 

The overall efficiency e is estimated from MC sim- 
ulations; it has a constant value of around 70% 
at energies > 3 TeV for all nuclei. 
The probability that a nucleus undergoes a spal- 
lation reaction in the amount of material (~ 4.8 
g/cm^) above the upper SCD plane is also es- 
timated from MC simulations. The fraction of 
surviving nuclei, i.e., the TOI correction, spans 
from 81.3% for C to 61.9% for Fe. The TOA cor- 
rection is calculated by simulating with FLUKA 
the atmospheric overburden during the flight (3.9 
g/cm^ on average). Survival probabilities ranging 
from 84.2% for C to 71.6% for Fe are found. 

5. Results 

The preliminary energy spectra of C, O, Ne, 
Mg, Si and Fe measured by CREAM-II are shown 
in figure [31 Only statistical errors are reported, 
a detailed study of the systematic uncertainties 
being underway. Absolute particle intensities are 
presented without any arbitrary normalization to 
previous data. The particle energy range extends 
from around 800 GeV up to 100 TeV. CREAM-II 
data are found in general to be in good agreement 
with measurements of previous balloon-borne and 
satellite experiments |7|8|9|10] . Though still pre- 
liminary, they seem to suggest that the intensities 
of the more abundant heavy elements have a very 
similar energy dependence. A more refined anal- 
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Figure 3. Energy spectra of the more abundant heavy nuclei. Results of CREAM-II (filled circles) are 
compared with measurements from HE AO [7 (triangles), CRN [8] (squares), ATIC [9 (open circles) and 
TRACER [IF (stars). 



ysis, including an assessment of the systematics, 
is still in progress. 

6. Conclusion 

The CREAM-II instrument carried out mea- 
surements of high-Z cosmic-ray nuclei with an ex- 
cellent charge resolution and a reliable determina- 
tion of their energy. Energy spectra of the more 
abundant heavy nuclei are measured and found 
to agree well with other direct measurements. 
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